Generalized priority-queue network dynamics: Impact of team and hierarchy 
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We study the eflect of team and hierarchy on the waiting-time dynamics of priority-queue net- 
works. To this end, we introduce generalized priority-queue network models incorporating interac- 
tion rules based on team-execution and hierarchy in decision making, respectively. It is numerically 
found that the waiting time distribution exhibits a power law for long waiting times in both cases, 
yet with different exponents depending on the team size and the position of queue nodes in the hi- 
erarchy, respectively. The observed power-law behaviors have in many cases a corresponding single 
or pairwise-interacting queue dynamics, suggesting that the pairwise interaction may constitute a 
major dynamic consequence in the priority-queue networks. It is also found that the reciprocity of 
influence is a relevant factor for the priority- queue network dynamics. 

PACS numbers: 89.75.Da, 02.50.Le, 89.65.Ef 



Introduction 

Priority-based queueing models are of interest among 
statistical physics community recently [ll, [2, H, 0, H, H, H 

. They were originally introduced in the context of hu- 
man dynamics [l|, aiming to model individual's internal 
decision making process underlying observed human ac- 
tivity patterns, yet bear implications to other problems in 
statistical physics such as the Bak-Sneppen model for bi- 
ological evolution ^ and invasion percolation [s!] as well, 
from the extremal dynamics perspective. Thus the un- 
derstanding of dynamics of the priority-queue models is 
of broader interest beyond the human dynamics. 

The primary fingerprint of priority-queue dynamics is 
the heavy tail in the waiting time distribution P(t), often 
taking an asymptotic power-law form for large waiting 
times. 



P(r) 



(1) 



where the waiting time r is the time interval between 
two consecutive events. In the original single fixed-length 
queue model by Barabasi [l| , the exponent a k. \ was ob- 
tained numerically, which was verified analytically later 
[^] . Such heavy-tailed dynamics have been observed for a 
range of human activities such as the e-mail, library loan, 
web browsing, and call initiation from a mobile phone 
[l, JH, XL, J^, JJi], prompting further interests and de- 
bates. 

Later on, the Barabasi model has been modified to 
include different factors such as memory and multi-task 
executions that would be of potential relevance in human 
decision making fT3. [Tsl. [igI. [l7j . Among them, one of the 
most confounding factor would be the human interaction 
[^,0]: In a modern society, the activity of an individ- 
ual is rarely an outcome of completely autonomous deci- 
sions, but of delicate compromises and balanced conflicts 
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between often competing priorities of various settings in 
complex social networks JJi, il9„ 2^ ■ Beyond its sociolog- 
ical context, the human interaction introduces coupling 
of multiple queue dynamics, rendering the model highly 
nontrivial. Thus the study of the role of human interac- 
tion on priority-queue dynamics is of theoretical interest 
as well. 

There are a few studies investigating the impact of in- 
teractions on the patterns of priority-queue dynamics. 
Oliveira and Vazquez introduced a model of two inter- 
acting priority queues with AND-type interaction 
and found that the waiting time distribution P{t) of the 
model is still a power law for large r, yet with different 
exponents such as a/ = 2 for the /-tasks in the L — 2 
case. Their model has been further generalized into the 
priority-queue networks of iV > 2 interacting queues 0] , 
where it is found that the AND-type interaction is not 
suitable for the network with loops when N > 2, leading 
to frozen, trivial queue dynamics. Therefore a scalable 
interaction of OR-type rule is introduced in the priority- 
queue network (the OR model Q). The OR model leads 
to different P{t) than the AND-type one and the power- 
law exponent a is found to depend on various factors 
such as the network size, global network topology, and 
local position of a queue node in a diverse way. These 
works have demonstrated that the interaction is indeed 
a relevant and consequential factor in the priority-queue 
dynamics. 

In this paper, we extend the work of Ref. Q, and con- 
sider further forms of human interaction in the priority- 
queue networks. Here we are specifically interested in two 
factors: the team-based task execution and the hierarchy 
in decision making. The former refers to the situation 
when a task demands simultaneous actions of more than 
two individuals and an individual's decision is affected 
by more than one others. The latter applies when there 
is a hierarchy in the queue nodes' status in a way that 
a node in higher hierarchy can order the execution of a 
task to the node in lower hierarchy. Both forms of human 
interaction are encountered in many real-life situations, 
thus the understanding of their impact is essential for 
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FIG. 1; (Color Online) Schematic illustration of the interac- 
tion rules of the models, (a) (A*", r)-model for A'' = 4 queues 
with teams of size r = 3. Each queue has (j) /-tasks (shaded) 
and one O-tasks (white). Each team is indicated by different 
shade level and labeled with team members, (b) ff-model for 
A'' = 4. Each queue has one /-task and one O-task. Arranged 
by the hierarchy value from top to bottom, when the a node, 
say the h = 2 node, executes its /-task, the nodes in the lower 
hierarchy, h — 3 and /i = 4 nodes, has to follow to execute 
their /-task (dark). 



of equal fixed length of Li = (^_{^) + 1 (Fig. la). In this 
model, when a node (chosen randomly in each step) has 
the J-task, Ig, as its highest priority task, all the nodes 
in the team G executes it at the step; an OR-type rule. 
The case with r = 2 is the same as the original OR model 
with the pairwise task execution of Ref. Q. As r > 2, 
the number of /-tasks &i changes drastically, which is ex- 
pected to affect the queue network's dynamics and thus 
is of interest in this work. This model can be thought of 
as the generalization of the priority-queue network model 
into the priority-queue hypergraphs. 

Interaction rule in the hierarchical decision — In the 
hierarchy-based task execution model (the //-model, 
hereafter) , N individuals are assigned its hierarchy value 
from h — 1 to h = N. Each node has two tasks, one 
/-task and an O-task; thus Li — 2 (Fig. lb). When a 
node i {i-th in the hierarchy) chooses its /-task to be the 
highest priority, all the lower hierarchy nodes (from i -I- 1 
to N) follow to execute it simultaneously with the node 
i. Depending on the node's position in the hierarchy, the 
degree of interruption (number of /-task calls from other 
nodes) varies, the effect of which is of interest in this 
work. 



a more complete human dynamics modeling, as well as 
for a more thorough understanding of the priority-queue 
network dynamics in general. 



Models 

The models studied in this paper are built upon the 
priority-queue network model of Ref. Q, with newly- 
introduced interaction rules. Depending on the interac- 
tion rule, each queue node has Qi /-tasks and an O-task 
{Li — Qi + 1). Initially each task is given a priority 
value drawn from a uniform distribution in [0, 1). Then 
each step, a queue node is chosen randomly (say it to be 
i) and its highest priority task is identified. If it is an 
/-task, its execution is challenged against the queue dis- 
cipline (interaction rules) of the particular model. If the 
selected /-task does pass the challenge, the node i and 
all other nodes involved in the task executes it simulta- 
neously. Otherwise, only the node i executes its O-task 
instead. Upon execution, the waiting time r of the tasks 
is recorded, and the executed tasks are replaced with new 
tasks each with a random priority value in uniform [0, 1). 
N such updates constitute a Monte Carlo step, which is 
the time unit of waiting time measurement. 

Interaction rule in the team-based task execution — In 
the team-based task execution model, each /-task is as- 
sociated with a group of queue nodes of size r, meaning 
that more than two individuals are involved in the ex- 
ecution of the task. Here the size of a team r is the 
important parameter characterizing the model. We call 
it an {N, r)-model hereafter. In the (A^, r)-model, each 
queue i has 9^ = -f-tasks plus an O-task, thus is 



Results 

Team-based task execution model: The {N,r)-model — 
In Fig. 2, we show the P{t) of the (A^, r)-model with 
N — 4,5,6. As mentioned above, r = 2 case of the 
{N, r)-model is identical to the OR model of Ref. 0] on 
the fully-connected network and we obtain that a/ ~ 2 
and ao ~ 1.3 (Fig. 2, □). For 3 < r < iV- 1, the waiting 
time dynamics is found to have the same asymptotics as 
the r = 2 case (Fig. 2, V, O, O). Notable exception is the 
/-task dynamics for r ~ N — I, for which the waiting time 
exponent is larger than the other cases as ai_r=N-i ~ 
2.5 (Figs. 2a-c, O). Interestingly, this distinct power-law 
exponent for r = A^ — 1 has also been observed in the 
pairwise (r = 2) OR model on fully-connected networks. 
There the case with N = 3 exhibits a distinct power- 
law exponent than other cases with A^ > 4 0]. The 
origin of this distinct behavior oi r — N ~ 1 case is not 
fully understood yet, calling for a further study. For the 
O-task dynamics, although the probability density for 
large r becomes elevated with r, the power-law exponent 
remains unchanged (Figs. 2d-f). 

The case with r = N exhibits a distinct behavior. Its 
/-task dynamics exhibits a fast decaying, exponential- 
like P(t) and the O-task dynamics shows a power-law 
with the exponent ao.r=N ~ 2, distinct from the cases 
with r < N (Fig. 2, A). For r — N, each queue node 
has one /-task and one O-task, similar to the N = 2 OR 
model. However, the /-task should be executed if any one 
queue has the /-task as its highest priority task, giving a 
disproportionately high probability to /-task executions, 
thereby leading to the exponential-like P{t) for the /- 
task. Meanwhile, from the perspective of the O-task, the 
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FIG. 2: (Color Online) Waiting time distribution P(r) of 
the team-based task execution model, the (A*', r)-model. (a— 
c) The waiting time distribution P(r) of /-tasks for various 
team size r with A'^ = 4 (a), 5 (b), and 6 (c). Slope of the 
straight lines is —2, drawn for the eye. (d— f) The waiting 
time distribution P{t) of 0-tasks for various team size r with 
= 4 (d), 5 (e), and 6 (f). Slopes of the straight lines are 
— 1.3 (upper) and —2 (lower), respectively, drawn for the eye. 



situation is not different from that of the N = 2 OR 
model, thereby leading to the same exponent ao ~ 2 as 
that of the TV = 2 OR model f?\. 

In the {N, r)-model, the mean /-task waiting time is fi- 
nite but the variance diverges for r < TV {aj < 3). For the 
O-task, not only the variance but also the mean waiting 
time diverge (a < 2), implying the lack of characteristic 
scale in the waiting times dynamics [2l'| . 

Hierarchical decision model: The H-model — In Fig. 3, 
we show the P{t) of the //-model with 2, 4, and 8. 
In the //-model, the node at the top of hierarchy {h — 1) 
is not affected by any other nodes in its decision making. 
Thus the h = 1 node's dynamics, having one /-task and 
one O-task, is the same as that of the Barabasi queue 
with length L = 2 Thus we have a = 1 for both the 
/-task and O-task for the h = 1 node (Fig. 3, □). 

The /-task dynamics of nodes with /i > 1 is found 
to follow the power-law P{t) with the same exponent 
a/ « 1, although these nodes' /-task dynamics is affected 
by the nodes in the higher hierarchy. This behavior can 




FIG. 3; (Color Onhne) Waiting time distribution P{r) of the 
hierarchical decision model, the /f-model. (a— c) The waiting 
time distribution P(r) of /-tasks for various hierarchy value 
h with TV = 2 (a), 4 (b), and 8 (c). Slope of the straight lines 
is —1, drawn for the eye. (d— f) The waiting time distribution 
P(r) of O-tasks for various hierarchy value h with TV = 4 (d), 
5 (e), and 6 (f). Slopes of the straight lines are —2 (upper) 
and —1 (lower), respectively, drawn for the eye. 



be understood by the following argument. In the H- 
model, a node's /-task dynamics is the superposition of 
subdynamics having power-law P(r) with tail exponent 
— 1, originating from the node itself and its higher hierar- 
chy nodes. Due to the singular nature of the power-law 
distribution with the exponent —1, the superposed dy- 
namics maintains the same tail exponent a ~ 1, upto the 
characteristic t which decreases with h, namely the num- 
ber of superposed subdynamics. Therefore, /-task dy- 
namics of the //-model follows the power-law P(t) with 
the same exponent a/ ~ 1 independent of node's position 
in the hierarchy. Meanwhile, although the power-law ex- 
ponent ai is close to 1, the mean waiting time is finite, 
as the mean value is dominated by the strong peak of 
P(r) at r = 1. However, the power-law tail with ai ~ 1 
means that once the /-task is shelved incidentally, it may 
wait for an extremely long time before execution. 

The O-task dynamics of nodes with /i > 1 is affected 
more significantly by the hierarchical decision. For these 
nodes, the relative priority of the O-task can be changed 
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by the forced execution of the /-task following the call 
from the higher hierarchy nodes. Thus the situation be- 
comes similar to the O-task dynamics in the OR model 
with N — 2 nodes having ao ~ 2, where the J-task pri- 
ority can be updated via the /-task call from the other 
node regardless of its own priority value. This similarity 
explains the exponent ao ~ 2 for the h > I nodes in 
the //-model. Whereas for the highest hierarchy {h = 1) 
node the mean O-task waiting time is finite, for h > 1 
nodes it is marginally diverging (ao ~ 2), suggesting a 
lack of the characteristic waiting time. 

Summary 

In this paper, we have considered two additional forms 
of interaction in the priority-queue network models, mo- 
tivated by the team-based task execution and hierarchy- 
based decision-making, which are encountered in many 
real-life human activities. Our numerical study has 
shown that these generalizations maintain the power- 
law decaying waiting time distribution P{t), Eq. (1), yet 
yield different values of the exponent a depending on the 
model parameters. 



Most of the power-law behaviors exhibited by the 
team-based {N, r)-model are found to have a counter- 
part in the pairwise OR model's dynamics. This suggests 
that the two-body interactions may contain the most 
fundamental dynamic consequence of interactions in the 
priority-queue network dynamics, although the precise 
mechanism of these correspondences has yet to be un- 
derstood, due to the lack of analytic results. For the 
hierarchical decision-based //-model, it is found some- 
what counterintuitively that it is the O-task dynamics, 
rather than the /-task dynamics, that is more strongly 
affected due to the hierarchical, directional interaction, 
compared to the reciprocal interaction case of the pair- 
wise OR model in which both the /- and O-task dynam- 
ics are affected. It thus highlights the nontrivial role of 
reciprocity of interaction in the priority-queue dynamics. 
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